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ABSTRACT 

Star formation is slow, in the sense that the gas consumption time is much longer than the dynamical time. It 
is also inefficient; essentially all star formation in local galaxies takes place in giant molecular clouds (GMCs), 
but the fraction of a GMC converted to stars is very small, ~ 5%. While there is some disagreement over 
the lifespan of GMCs, there is a consensus that it is no more than a few cloud dynamical times. In the most 
luminous starbursts, the GMC lifetime is shorter than the main sequence lifetime of even the most massive 
stars, so that supernovae can play no role in GMC disruption; another feedback mechanism must dominate. 
We investigate the disruption of GMCs across a wide range of galaxies, from normal spirals to the densest 
starbursts; we take into account the effects of HII gas pressure, shocked stellar winds, protostellar jets, and 
radiation pressure produced by the absorption and scattering of starlight on dust grains. In the Milky Way, 
we find that a combination of three mechanisms — jets, HII gas pressure, and radiation pressure — disrupts 
the clouds. In more rapidly star forming galaxies such as "clump" galaxies at high-redshift, ultra-luminous 
infrared galaxies (ULIRGs) and submillimeter galaxies, radiation pressure dominates natal cloud distribution. 
We predict the presence of 10-20 clusters with masses ~ IO^Mq in local ULIRGs such as Arp 220 and a 
similar number of clusters with M» ~ IO^Mq in high redshift clump galaxies; submillimeter galaxies will have 
even more massive clusters. We find that the mass fraction of a GMC that ends up in stars is an increasing 
function of the gas surface density of a galaxy, reaching ^ 35% in the most luminous starbursts. Furthermore, 
the disruption of bubbles by radiation pressure stirs the interstellar medium to velocities of ~ 10 km s~' in 
normal galaxies and to 100km s~' in ULIRGs like Arp 220, consistent with observations. Thus, radiation 
pressure may play a dominant role in the ISM of star-forming galaxies. 



Subject headings: Galaxies: star clusters, formation, general, starburst - 
— stars: formation 



HII regions — ISM: clouds, bubbles 



1. INTRODUCTION 



The Kennicutt law (Kennicutt 1998) 



(1) 



relates the star formation surface density E* to the gas sur- 
face density and the local dynamical time ft w Vc/R in 
disk galaxies, where Vc is the circular velocity of the galaxy 
and R is the distance from the galactic center The dimension- 
less constant rj w 0.017 is surprisingly small, a finding that 
is interpreted as showing that star formation is a slow pro- 
cess. Star formation is similarly slow on smaller scales within 
galaxies (Kennicutt et al. 2007; Bigiel et al. 2008; Leroy et 
al. 2008; Krumholz & Tan 2007). 

Remarkably, equation ([T]l holds for galaxies like the Milky 
Way, with rather modest star formation rates of order a so- 
lar mass per year, for starburst galaxies with star formation 
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rates of order tens of solar masses per year, for ultraluminous 
infrared galaxies (ULIRGs) with star formation rates around 
one hundred solar masses per year, and for sub-millimeter 
galaxies with star formation rates in excess of one thousand 
solar masses per year. There are indications, however, that 
the star formation efficiency 77 may be larger in UL IRGs and 
sub-mm galaxies, with 77 0.1 (iBouche et al.ll2007h . 

The large range of galaxies that obey equation ^ suggests 
that whatever process sets the efficiency of star formation op- 
erates in galaxies with very different conditions in their in- 
terstellar media. For example, the gas surfa ce density in the 
Milky Way at 8kpc is « 2 x 10"''gcm"2 ( Boulares & Coxl 
II99OI) . while that in the two lOOpc star forming disks of the 
ultra-luminous infrared galaxy (ULIRG) Arp 220 is T,g w 
7 gem"-; the mean gas densities of the two galaxy's star form- 
ing disks also differ by a factor of > 10"^. Although the range 
in turbulent velocities in the ISM is not so dramatic, from 
6km s"' in the Milky W ay to 60- 80km s"' in Arp 220 



jDownes & Solomonll998 ). the turbulent pressure in Arp 220 
exceeds that in the Galaxy by a factor of ~ 10^. 

Another signature of inefficient star formation relates to 
individual giant molecular clouds (GMCs). In the Milky 
Way, all stars are believed to form in such clouds. How- 
ever, the fraction egmc of a GMC that is turned into stars is 
quite l ow, around 5% in the Milky Way (IWilliams & McKe3 
,1997i; lEvans et al.ll2008h . with a similar value inferred from 
the global star formation efficiencies in other galaxies (see. 
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e.g. jMcKee & Ostrikerfl2007h . 

One class of explanation for this low star formation effi- 
ciency is that gas in the ISM is prevented from collapsing by, 
for example, magnetic fields, cosmic rays, or by externally- 
driven turbulence (Parker 1969; Sellwood & Balbus 1999; 
Ostriker, Stone, & Gammie 2001). A second class of explana- 
tion is known by the name of "feedback:" the injection of en- 
ergy and momentum into the ISM by stellar processes so that 
star formation alters the ISM conditions and limits the rate at 
which gas turns into stars. The form the feedback takes is not 
currently agreed upon. Suggested mechanisms include super- 
nova heating, deposition of momentum by supernovae, heat- 
ing by photoionizing radiation from massive stars, deposition 
of momentum by expanding bubbles of photoionized HII re- 
gion gas, deposition of momentum by the shocked winds from 
massive stars, and jets from protostars (e.g., McKee & Os- 
triker 1977; Silk 1997; Wada & Norman 2001; Matzner 2002; 
Li & Nakamura 2006; Cunningham et al. 2008). 

In this paper, we study these feedback processes and assess 
the role that they play in disrupting GMCs across a wide range 
of star-forming galaxies. In addition, we focus on a some- 
what less well studied form of feedback: deposition of mo- 
mentur n by the absorption and scattering of starhght by dust 
grains dOlD ell et al.' 1967; 'Chiao & Wickramasinghe 1972'; 
Elmegree"nl,1 983; Ferrara 1993; Scoville et al. 2001; Scoville 
2003UThompson et all 120051) . Although the magnetic fields 
in starburst galaxies can be large (^few mG for Arp 220; 
Thompson et al. 2006, Robishaw et al. 2008), we neglect them 
throughout this paper in order to focus on the competition be- 
tween various processes that contribute to disrupting GMCs. 

It is important to distinguish between two arenas in which 
galactic feedback likely operates: galactic disks in the large, 
and in the main units of star formation, GMCs. While there 
is rather sharp debate in the literature, we will assume that 
GMCs are at least marginally gravitationally bound objects, 
and hence that they are unlikely to be supported by feedback 
acting on the scale of galactic disks as a whole. As noted 
above, observations in our own and nearby galaxies estab- 
lish that only ~ 5% of the gas in a GMC ends up in stars; 
the rest of the gas is dispersed back into the ISM. The uni- 
versality of the Kennicutt law suggests that a similarly small 
fraction of the GMCs in other, more distant classes of galax- 
ies, is turned into stars. Something is disrupting GMCs, but 
it is unlikely to be large scale turbulence in the galaxy as a 
whole. Instead, GMCs must be disrupted by the stars that 
form in them. A number of authors have argued that galactic 
GMCs are disrupted by expanding HII regions (e.g., Matzned 
12002; Krumholz et al. 2006); this mechanism cannot, how- 
ever, work in luminous starbursts (Matzner 2002). The fact 
that these galaxies nonetheless have roughly similar star for- 
mation efficiencies suggests that another disruption mecha- 
nism must be competitive with expanding HII regions. 

In this paper, we argue that the radiation pressure produced 
by the largest few star clusters in a GMC, acting on dust grains 
in the gas, is the primary mechanism by which GMCs are 
disrupted in more luminous starbursts and massive GMCs (see 
also Scoville et al. 2001, Harper-Clark & MuiTay 2009, and 
Pellegrini et al. 2007, 2009). Protostellar jets also provide 
an important contribution, particularly in the early stages of 
the evolution. In spirals like the Milky Way, both expanding 
HII regions and radiation pressure are comparably important, 
depending on the size and mass of the cluster, and supernovae 



also play an important role in the latest stages of disruption. 

1.1. Is Feedback Really Necessary? 

A key thesis of this work is that stellar feedback is crucial 
for understanding the low observed val ues of the star forma- 
tion ef ficiency in galaxies. In contrast, iKrumholz & McKe3 
(l2005h present an explanation of the Kennicutt-Schmidt law 
(eq. dJ that does not invoke an explicit form of feedback. 
Their argument is that turbulent motions prevent the collapse 
of the bulk of the gas in a GMC (or in other bound objects); 
only if the density is above a critical density, which depends 
on the Mach number of the flow, do stars actually form. The 
fraction of gas in a turbulent flow that lies above this critical 
density is small, leading to the low observed star formation 
efficiency per dynamical time. 

We find this argument to be compelling, as far as it goes. 
As long as turbulence is maintained, only a small fraction 
of gas will collapse into stars per dynamical time. How- 
ever, the ass umption of a constant level o f turbulence is es- 
sential to the IKrumholz & McKe3 (l2005h argument. A key, 
and yet unanswered, question is thus what maintains the tur- 
bulence? If the turbulence in a GMC is not maintained, then 
the GMC will contract, leading to an increase in the mean 
density and a decrease in the dynamical time. Indeed, sim- 
ulations find tha t turbulence decays on ~ 1 crossing time 
jMac Low! 119991; lOstriker et"aLl 1200 lb . so fliat a continued 
source of energy is needed to maintain the turbulent support 
of the cloud. It is possible, in principle, that gravitational 
contraction of a GMC can sustain the turbulence, maintain- 
ing appro ximate virial equ ilibrium and a slow contraction of 
the cloud (IKrumholz et al. 2006). We argue, however, that an 
independent internal source of turbulence, provided by stars, 
is crucial to maintaining the slow rate of star formation. 

As an example, we apply this argument to Arp 220. The 
interstellar medium of Arp 220 has a turbulent Mach number 
~ 100. The fraction of a GMC (or any bit of molecular 
gas) that is sufficiently dense to be converted into stars in a 
free fall time is then ~ 0.013-0.05 for GMCs with a virial 
param eter = 0.1-1 (see Figure 3 of Krumholz & McKe^ 
l2005h .^ According to this argument, a GMC will convert half 
its gas into stars in ten to forty free fall times, reasonably 
consistent with equation ([T]| for any avii. However, this as- 
sumes that the cloud does not collapse and reduce its free fall 
time. In reality, if turbulence can only maintain avir '--^ 0.1, 
a cloud is likely to collapse, leading to a rapid increase in 
density and a decrease in the free-fall time. If the star for- 
mation rate per free-fall time remains roughly constant, then 
the actual star formation rate will increase rapidly with time, 
and most of the gas in the cloud will be turned into stars in 
roughly on e initial (la rge-scale) free-fall time. Thus the model 
of Krumholz & McKee ( 2005) for the low star formation ef- 
ficiency in galaxies relies critically on maintaining sufficient 
levels of turbulence so that avir ^ 1 . On larger scales — above 
the characteristic GMC size — the equivalent argument is that 
the galactic disk must have Toomre 2 1, as we discuss be- 
low. 

1.2. This Paper 

^ A viral pai'ameter ayii- = 1 corresponds to a cloud that is just gravitation- 
ally bound, while the smaller value of Ovir = 0. 1 corresponds to a cloud that 
is approaching free fall conditions. 
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The remainder of this paper is organized as follows. In ^ 
we collect a number of relevant astrophysical results used in 
our modeling. In ^ we describe a simple one-dimensional 
model for the disruption of GMCs which includes the ef- 
fects of HII gas pressure, proto-stellar jets, radiation pressure, 
gas pressure associated with shocked stellar winds, and wind 
shock generated cosmic rays (many of the details of how we 
model these forces are given in Appendix IaIi. In section ^ 
we present the results of our numerical modeling of GMC 
disruption in star forming galaxies. To explore the wide range 
of conditions seen in galaxies across the Kennicutt-Schmidt 
law, we consider models for GMCs in the Milky Way, M82, 
Arp 220, and the z - 2 galaxy Q2346-BX 482. In ® we dis- 
cuss the implications of our results, the origin of turbulence in 
galaxies, and explain physically why radiation pressure is the 
only source of stellar feedback in principle capable of disrupt- 
ing GMCs across the huge dynamic range in ISM conditions 
from normal galaxies to the densest starbursts. 

2. ASTROPHYSICAL ELEMENTS 

In this section we collect several pieces of observations and 
physics that we believe are relevant to star formation in galac- 
tic disks and GMCs. We order these items according to the 
amount of support they enjoy, from substantial to slim. The 
key conclusion below is that a significant fraction of all stars 
are formed in compact {^^ few pc radius) massive star clusters, 
that in turn reside in GMCs. Given the importance of a few 
star clusters that are small compared to the GMC as a whole, 
a one dimensional model for GMC disruption is a reasonable 
first approximation; this is presented in ^ 

2.1. Marginally Stable Disks (Q~ I) 

Quirk (1972) showed that normal galaxies have gas disks 
with 2 ~ 1- Kennicutt (1989) refined this to the statement 
that within the star forming part of normal galactic disks, 
1/4^2^ 0.6. At large radii he found Q > 1 and a lack 
of star formation. More recently, Leroy et al. (2008) studied 
star formation in detail in 23 nearby galaxies. They found 
that if they accounted for only the gas surface density, as done 
above, their disks were stable, with g « 3-4; using the total 
(gas plus stars) surface dens ity resulted inj2_~ 2 with a slight 
variation in radius. Unlike iKennicuta (Il989l) . they find that 
star formation occurs at large radii, beyond Kennicutt's Q> I 
radius, albeit at reduced rates. 

These studies were restricted to normal galaxies, and em- 
ployed a fixed sound speed as the estimate for the random 
velocity. However, there is evidence that disks with turbulent 
velocity vt > > Cj also satisfy Q ^ 1 when vj- is used in evalu- 
ating Q (e.g.. [Thompson et al.l2005 and our discussion of the 
starbursts M82, Arp 220, and Q2346-BX 482 in ® . Moti- 
vated by these observations and by theoretical considerations, 
we will assume that all star forming galaxies have 2 « 1 . 



This scenario is consistent with observations of GMCs in 
our galaxy; in the Milky Way half the gas is in molecular 
form in giant mo lecular clouds with a characteristic mass of 
order 5 x 10^ Mq (ISolomon et al.l 19871) . but with a rather wide 
range of masses. The number of clouds N(m) of mass m is 
given by 

f^=A.o(^)""^ (2) 

am \ m / 

with an exponent ag ^ 1.8 (ISolomon et alj [T987h or 1.6 
dWilhams & McKe^ Il997h . so that most of the mass is in 
the largest clouds. A s a cautionary note, we note that 
lEngargiola eTal] (l2003h find ac = 2.6 ± 0.3 in M33, which 
suggests that lower mass clouds contribute a significant frac- 
tion of the total mass. In the Milky Way the largest GMCs 
have masses of order - 3 x lO^'Mg (ISolomon et all [1987). 
roughly consistent with the Toomre mass. In the Milky Way, 
and possibly in other galaxies, molecular clouds are sur- 
rounded by atomic gas with a similar or slightly smaller mass. 

The clouds appear to be somewhat centrally concentrated. 
We will often employ a Larson-law density distribution. 



p(r)' 



(3) 



where r is the distance from the center of the GMC. We also 
explored isothermal models p(r) 1 /r^; we find that such 
clouds are slightly easier to disrupt than the less centrally con- 
centrated Larson-law clouds in the optically thick limit. 

2.2.1. But are there molecular clouds in ULIRGs? 

GMCs are observed in the Milky Way, and in nearby star- 
forming galaxies such as M82. We see clumps of gas in 
"chain" or "clump" star-forming galaxies at z = 2, such as 
Q2346-BX 482 discussed below. These have been inte rpreted 
as self-gravitating, i.e., as GMCs (e.g.. iGenzel et alj |2008). 
However, we do not know of any direct evidence for Toomre 
mass self-gravitating objects in ULIRGs. There is some evi- 
dence against such objects: since the clouds are self gravitat- 
ing, they will have a slightly higher velocity dispersion than 
that of the disk out of which they form. Increasing the ve- 
locity dispersion will alter the inferred gas mass (see Downes 
& Solomon 1998). Putting too much gas in gravitationally 
bound objects will increase the apparent gas mass, possibly 
making it larger than the dynamical mass. 

One the other hand, ongoing star formation is clearly seen 
in Arp 220. Star formation probably requires densities ex- 
ceeding 10'' cm"^^ to proceed. Thus there is evidence that 
some gas is gravitationally bound. Moreover, there are nu- 
merous massive co mpact star clusters observed in Arp 220 
dWilson et alj|2006h . indicating that massive, bound, and rel- 
atively compact accumulations of gas existed in the recent 
past. Motivated by these considerations, we will assume that 
Toomre-mass self-gravitating objects exist in all star-forming 
galaxies, including ULIRGs. 



2.2. The Toomre Mass and Giant Molecular Clouds 

We assume that galactic disks initially fragment on the disk 
scale height H « (vy / Vc)r. The fragments will form gravi- 
tationally bound structures with a mass given by the Toomre 
mass, Mt — nH^J^g. Near the location of the sun, the gas 
surface density « 2 x 10"^ gem"- and H ~ 300pc, giving 
Mt~2x IQ^Mq. 



2.3. Gas clump and stellar cluster mass distributions 

Most of the gas in Milky Way GMCs is diffuse (n < 
3 X 10^ cm"^), but a fraction of order 10 % is in the form of 
dense gas clumps, with sizes around 1 pc dLada & Ladal2003h 
and masses from a few tens to a few thousand solar masses. 
The clumps have a mass distributio n similar to that for clouds 
(eq.|2l), with an exponent ac ~ 1.7 dLada & Ladall2003l) . 



4 



TABLE 1 

OBSERVED GALAXY PARAMETERS 



Galaxy 




^dyn ^ 




Z/Zq 


M, (obs) 




kpc 


yr g cm"^ 


km S-' 




Mq yr-' 


Milky Way 


8.0 


3.6 X 10^ 2 X 10-3 


6 


1 


2.0 


M82 


0.35 


3.0x106 Qi 


15 


1.5 


4 


BX482 


7.0 


2.9 X 10' 4 X 10-2 


53 


1 


140 


Arp220 


0.1 


3.3 X 10^ 7 


61 


3 


120 


Note. — 


Observed galaxy properties. Column one 


gives the name of the 



model. The next five columns give model input parameters: the disk radius 
(col. 2), dynamical time Rd/vc (col. 3), gas surface density (col. 4), turbulent 
velocity vt (col. 5; recall that H = [vt lvc]Rd), and metallicity in solar units 
(col. 6). The metallicity is not that well-constrained in Arp 220. Column 7 
gives the observed star formation rate for the galaxies. 



TABLE 2 

GMC AND STAR CLUSTER PROPERTIES 



Galaxy 


^GMC 


Mguc 




M, 


<:gmc 


Vt 




pc 


Mq 


pc 


Mq 




km S-' 


Milky Way 


100 


3 X 10* 


2 


10= 


0.03 


5 


M82 


23 


3 X 10*5 


1.5 


7 x 10' 


0.24 


10 


BX482 


925 


10' 


13 


2.7 x 10* 


0.27 


50 


Arp 220 


5 


4x 10' 


3.5 


1.4 X 10' 


0.38 


50 



Note. — Columns 2-5 give the assumed GMC and star cluster prop- 
erties: the radius of the GMC Rgmc (col. 2), the mass of the GMC 
Mqmc (col. 3), the star cluster radius (col. 4), and the stellar mass 
of the star cluster M, (col. 5). Columns 6 and 7 give the predictions of 
our model for the star formation efficiency in the GMC eoMC and the 
shell velocity when the GMC is disrupted, which we also interpret as 
the turbulent velocity vj induced in the ISM of the Galaxy. 



In both the Milky Way (Hlmegreen & Efremovl 119971: 
Ivan den Bergh & Lafontaine' '19841 ) and in nea rby galaxies 
iMcKee & Williams 1997; Kennicu ttetanil989h . the number 
of stellar clusters of mass m» is given by 



*cl 



*cl.O 



'«*cl 



(4) 



with ttci ~ 1.8. In other words, most stars form in massive 
clusters; in the MiUcy Way, at least, these clusters are made 
from gas in massive gas clumps, inside of massive GMCs. 

2.4. The Sizes of Star Clusters 

Star clusters are observed to have sizes ranging from rd ~ 
0.1 pc (forMd « IOMq) to rd « lOpc (Md « IO^Mq). There 
are hints that clusters with masses Md < lO^^M© have a mass- 
radius relation of the form 



Md\ 
mo / 



/3 



pc 



(5) 



with mo = IO'^Mq and (3 « 0.4 dLada & Ladal l2003h . but 
it is entirely possible that this is a selection effect. In- 
termediate mass clusters, those with lO^M© < Md ^ 3 x 
lO^M© have rd ~ 2 pc independent of mass {(3 = 0), al- 
beit with substantial scatter This characteristic size is seen 
for young (< 30Myr) and old (> 30Myr) clusters i n M51 
with masses in the range lO^-lO^M© (Scheepma ker et alj 
l2007h . f or super star clusters i n M82 with Mh = 10 ^-4 x 
IO^Mq (iMcCradv et all I200I iMcCradv & GrahamI 
and in globular clusters with masses ~ 10^- lO^M© 
[T996h . Finally, high -mass clusters with Md > IO^Mr have 
/3 = 0.6 and mo = lO^Mtp in equation ^ ( Walcher et al. 2005; 
Evstigneeva et al. 2007t iBarmbv et ail [20071: IReikuba et al.l 
2007; Murray 2009). When using the radii of stellar clus- 



ters in our GMC models, we will be guided by these observed 
mass-radius relations. 

3. A MODEL OF CLUSTER & GMC DISRUPTION 

We start by specifying the properties of the disk in which 
the GMC lives; the effective disk radius R^, the circular ve- 
locity Vc, the disk scale height Hd, the disk gas mass Md and 
metallicity Z/Zq (relative to solar). We choose the values of 
these parameters to match those of four galaxies to which we 
compare our models: the Milky Way, M82, Q2346-BX 482, 



and Arp 220. In that sense these are not free parameters. Table 
1 summarizes the observed input parameters for the galactic 
disks in the systems we model, while Table 2 gives the in- 
ferred or assumed input parameters related to the GMC and 
its central star cluster: the mass Mgmc and radius /?gmc of the 
GMCs and the stellar mass and radius of the star clusters. Ta- 
ble 2 also lists the GMCs star formation efficiency ecivic and 
the shell velocity when the GMC is disrupted - we interpret 
the latter as the turbulent velocity induced in the ISM, vj- 

We employ a one-dimensional model for the GMC, with 
the free parameter 0g = Hd /Rguc defining the size of a GMC 
with respect to the disk in which it resides. In the Milky Way 
this ratio is about ^2-5 for the most massive clouds; for 
specific cases like G298.4-0.3, which has /?gmc ~ lOOpc, we 
use the observed ratio. In models for other galaxies we fix 
4>G = 4. The mass of the GMC is taken to be the Toomre 
mass, with a Larson-like (p ^ 1 /r) internal density profile. 

For the purposes of our simplified modeling, the stars are 
assumed to lie in a single massive cluster of total mass (gas 
plus stars) Md, which forms a mass of stars M* = edMd with 
luminosity L, surrounded by the remnants of the gas out of 
which it forrned, vv ith mass Mg = (1 -CciWci- We use a 
iMuench et alj (|2002|) stellar initial mass function (IMF) to re- 
late the cluster luminosity to its mass. The quantity ed char- 
acterizes the efficiency with whic h cluster gas is turn ed into 
stars. In our galaxy eri 0.3-0.5 dLada & Ladal2003l) . In our 
models, we fix EcI — 0.5 and adjust the cluster mass Md (or 
equivalently, the cluster stellar mass M*) to find under what 
conditions the central star cluster can disrupt its host GMC. 
Physically, we expect that star formation will self-adjust to a 
form a cluster of approximately this mass. For the Milky Way 
and M82 the star cluster masses we infer by this method are 
comparable to those observed. 

We model the impact of the central star cluster on the sur- 
rounding GMC using the thin shell approximation. As the 
star cluster evolves, driving winds, jets, and radiation into the 
overlying shell of gas, we calculate the shell's dynamics as it 
sweeps up mass and disrupts the GMC. The relation between 
shell radius r, shell velocity v(r), shell mass M(r), and shell 
momentum Psh('") is given by 



dr 
dt 



M(r)' 



(6) 



where M(r) is the mass of the shell, which increases as the 
shell radius r increases, and Psh = M(r)v(r) is the momentum 
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of the shell. For a Larson-Uke density profile of the GMC, 
p oc r~' , the mass of the shell is given by 



M(r) = (l-eci)Md+MGMC 



for r < Rgmc and by 



2 / \ 2 



Rgmc J V^gmc 



4tt 



M(r) = ( 1 - eci)Mci +Mgmc + y (''^ " ''^GMc)Pdisk 



(7) 



(8) 



for RcMC < r < H. We have experimented with both isother- 
mal (Mgmc('") ^ r) and Larson (Mgmc('") ^ r^) GMCs and 
find qualitatively similar results. 

The momentum equation for the shell is 

= --Pgrav - ^i-am " -Ptuib + -PhII + -Prad + -Pjet + -Fhot + -Fcr (9) 

The inward forces on the right-hand side of equation ^ are 
the self gravity /Vjav of the swept up shell and the mutual grav- 
ity of the stars in the largest cluster and the shell, the force Ftuib 
exerted by the turbulent motions of the gas in the GMC on the 
shell, and the force /Van, produced by ram pressure as the shell 
sweeps up the material of the GMC or the surrounding gas 
disk. The outward forces on the right-hand side of equation 
(|9]l acting to disrupt the GMC include the force /^et associated 
with momentum deposited by jets from star formation, the 
force Fmi due to the thermal pressure from ionized gas (HII 
regions) and from shocked stellar winds, and the force /vaj as- 
sociated with radiation pressure. In some models we include 
the effects of hot gas fhot and cosmic rays associated with 
shocked stellar wind. A more detailed description of how we 
implement these forces is given in AppendixlAl 

We assume that stellar wind energized hot gas lO^K) 
does not play a dynamical role in galaxies like the Milky 
Way; we show explicitly in this paper that even in the most 
optimistic cases such winds are not important for GMC 
disruption in ULIRGs. Much of the literature on bubbles 
around massive stars and massive star clusters assumes 
that shocke d stella r winds domi n ate th e dynamics (e.g.. 
Castor et al.1 119751: IWeaver et all 119771: IChu&Mac Lowl 
1990^ . However, observations of HII regions in the Milky 
Way suggest that the pressure in such hot gas is equal to that 
of th e associ ated HII (lO'^K) gas (Do rl a n d & M ontmerle 

T987t iMcKee. Van Buren. & LazareffI 119841: 

Harper-Clark & Murray! l2009l) . The most hkely inter- 
pretation of these results is that neither hot gas nor cosmic 
rays are confined inside bubbles in the Milky Way or the 
LMC, but rather escape (Harper-Clark & Murray 2009|). 
Accordingly, we neglect the last two terms on the right hand 
side of equation (|9]l in our Milky Way models and in our M82 
models; calculations which include these pressures predict 
bubble sizes that are far too large compared to observations 
in the Milky Way. Accordingly, with the exception of ^4.41 
we ignore the pressure associated with shocked stellar winds. 
We do include stellar-wind and cosmic -ray pressure in our 
Arp 220 models, but there they make little difference. 

In ^ we made a distinction between feedback in the disk 
as a whole, and feedback in GMCs. This distinction is im- 
portant for a number of reasons. In particular, we believe that 
supernova (SN) explosions largely contribute to the former, 
but not the latter In z ULIRGS, where the bulk of the 
star formation takes place in ~ lOOpc disks, the dynamical 



time is fdyn = R/vc ~5 x 10^ yr, much less than the main se- 
quence lifetime of even the most massive stars. Hence a grav- 
itationally bound object in a ULIRG cannot be disrupted on 
a dynamical time by SNe resulting from stars formed in that 
object. In the Milky Way, even though the dynamical time 
is longer than the main sequence lifetime of massive stars, it 
is clear that GMCs are in the process of bein g disrupted well 
before SNe occur (iMurrav & Rahmanll2009l) . Although SNe 
may be important during late stages of GMC disruption, and 
for stirring up the galactic disk as a whole, they cannot be 
the main agent that disrupts GMCs, either in Milky Way-like 
galaxies or in ULIRGs. For this reason, we do not include the 
force due to SNe in equation (|9]l. 

In the densest starbursts, which have mean gas densities 
~ 10''"'^ cm"^, SNe rapidly lose the majority of their energy 
to radiative losses (e.g., Thornton et al. 1998). Under these 
conditions, the primary role of SNe is to stir up the bulk of 
the ISM via the momentum they supply, rather than to heat 
up the gas and/or create a hot phase of the ISM in pr essure 
equilibrium with the rest of the mass ([Thompson et al.i r2005). 
To see that the latter is untenable, we estimate the density 
«/, of hot gas required for a virialized ISM at ^ lO^Tj K to 
be in pressure equilibrium with the bulk of the gas, i.e., for 
Ph = iihkT ~ vrGEg. Given this density, we find that there 
is a critical surface density above which thermal X-ray 
emission from the hot gas would exceed the observed cor- 
relation between X-ray and FIR emission {Lx — 10~^Lfir; 
iRanalh et al.ll20()3h : ~QM{Hd/lQQ^c)-^I^T^'\ Galaxies 
with Ed > Ec, which includes the majority of luminous star- 
forming galaxies, cannot have a dynamically important hot 
ISM. Instead, the dominant role of SNe is to stir up the dense 
gas via the momentum imparted in the snowplow phase. This 
may even be true in galaxies with E^ < E^/, because the hot 
gas can vent via galactic winds or fountains. For example, 
in the Milky Way, which has E,, ~ 2.5 x 10"^ g cm"- < E^, 
the hot ISM is believed to contribute only ^ 10% of the to- 
tal pressure (Boulares & Cox 1990). For this reason, we will 
not consider the pressure due to the hot ISM in this paper, 
although the momentum supplied by SNe is important in the 
late stages of star cluster and GMC evolution. 

Having summarized the basic elements of our model, we 
now describe its application to galaxies ranging from the 
Milky Way to the most luminous starbursts (®. We then 
discuss the implications of these results (©. 

4. RESULTS 

4.L The Milky Way: G298.4-0.3 

Recent work has revealed that the Milky Way harbors a 
number of very young massive star clusters, with rnasses 
ranging up to lO^M^ CFig eret al.ll l999HBrandner e"tan i2008l: 
iFiger et al.ll2006t iMurrav & Rahma^ |2009|) . We examine the 
case of G298.4-0.3 in the Carina arm. [Murray & RahmanI 
(2009) show that the free-free flux emerging from this region 
implies a total ionizing flux Q = 1.1 x 10^' s"', and suggest 
that ~ 60% of this comes from a single massive cluster re- 
siding in the prominent ^ 50pc bubble revealed by Spitzer 
GLIMPSE images. Most of the remaining flux comes from 
two clusters associated with the two giant HII regions G298.2- 
0.3 and G298. 9-0.4; both sources appear to lie in the rim of the 
bubble seen in the GLIMPSE images. Grabelsky et al.1 (1 19881) 
find two massive GMCs in this region, numbers 24 and 26, 
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Fig. 1 . — Shell radius as a function of time and velocity, forces, and shell momentum as a function of radius in our model for G298.4-0.3 in the Milky Way. 
The shell is sheared apart when it reaches the Hill radius (~ 200 pc), where we end our integration. Upper left: The dotted, dashed, and solid lines mark when the 
first cluster SNe explode, the central cluster luminosity drops to 1/3 of its initial value, and t reaches the Milky Way dynamical timescale, respectively. Upper 
right: The velocity of the swept-up shell in G298.4-0.3 as a function of shell radius. Note that the asymptotic velocity is comparable to the turbulent velocity of 
the Milky Way disk. Lower left: The upper-most solid line is the total outward force, consisting of the momentum supplied by protostellar jets (solid line; i^ei). 
by HIl gas pressure (short dashed line; fail) and radiation pressure on dust grains (long dash line; Frad). The dotted line is the total inward force (fi„), dominated 
by the self-gravity of the shell. Lower right: Momentum of the swept-up shell in G298.4-0.3 (solid line), together with the momentum deposited by radiation 
(long dashed line), gas pressure (short dashed line), and protostellar jets (solid line, labeled Pjei). The bulk of the momentum is supplied by radiation and gas 
pressure, but the early contribution of the protostellar jets is important in the disruption of the natal cluster gas. 



both with Mgmc « 3 x IQ^Mq and Rquc ~ lOOpc. Their ra- 
dial velocities are 22km s~' and 24km s~', in good agreement 
with the range of radio recombination line, i.e., HII region, ra- 
dial velocities in this direction, which range from -l-16km s~' 
to 30.3km s~', with a mean ~ -1-23 km s~'. 

Accordingly, our model for G298.4-0.3 consists of a GMC 
with /?GMC = lOOpc and Mgmc = 3 x IQ^Mq. In the spirit of 
our simplified one-dimensional modeling, we lump all of the 
star clusters together into a central star cluster with L « 7 x 
IO^Lq and initial cluster radius = 1.5pc. Half the Galactic 
star formation takes place in 17 star clusters, with a minimum 
2 = 3 X 10^^ so G298.4-0.3 is representative of star forming 
clusters in the Milky Way. 



The top two panels of Figure [T] show the radius of the shell 
surrounding the central cluster as a function of time and the 
shell velocity as a function of radius. The shell starts at our 
putative initial radius for the cluster of ^ 1.5pc, and reaches 
r ^ 80pc at about 6.5 Myr (dashed line), at which point the 
star cluster luminosity has dropped by a factor of 3. The most 
massive stars begin to explode after about 3.6 Myr (vertical 
dotted line), while the last O stars explode after about 1.3 x 
10^ yr. The solid vertical line marks the dynamical time R/vc 
for the Milky Way at 7? = 8 kpc. Note that the dynamical time 
for the GMC is somewhat shorter, ^ 6 Myr 

The lower left panel of Figure [T] shows the forces as a func- 
tion of radius in this model; note that the radiation and HII gas 
pressures drop after 6.5 Myrs, when the bubble has r ^ 100 
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pc, but the bubble continues to expand at the same rate. The 
evolution after 6.5 Myrs, i.e., radii larger than 50pc, may 
underestimate the rate of expansion somewhat; the bubble 
may expand slightly more rapidly after several Myrs due to 
energy input by SNe. We say the rate of expansion may be 
underestimated since the hot gas from the SNe is likely to es- 
cape the bubble as easily as the hot gas from shocked winds 
apparently does. In addition to SNe, other unmodeled effects 
also become important at late times and large radii. For ex- 
ample, because the inner radius of the shell exceeds the outer 
radius of the initial GMC, the surface density of the gas de- 
creases to Av w 1, so that ionizing photons from the Galactic 
radiation field can penetrate and ionize the shell. 

We halt the integration in our model when the radius of 
the expanding bubble exceeds the Hill radius rniii of the 
GMC, i.e., when the tidal shear from the Galaxy exceeds the 
self-gravity of the GMC: rniii ~ (MGMc/2M(r))'/''fl, where 
M(r) = v^a /G is the enclosed dynamical mass of the galaxy 
at the galactocentric radius a of the GMC. After this time the 
remaining molecular gas will be dispersed (although not nec- 
essarily converted to atomic gas). 

Figure[T]shows that the central cluster in G298.4-0.3 should 
disrupt its natal GMC. What force is responsible for this dis- 
ruption? At the current radius, ^ 55 pc, the radiation force 
and the gas pressure force are within a factor of two of each 
other, and will remain so until most of the O stars explode; the 
force from protostellar jets is substantially smaller However, 
at early times, the jet force was as much as a factor of two 
larger than the radiation pressure force, and the gas pressure 
force was negligible. 

Finally, the lower right panel of Figure [T] plots the mo- 
mentum of the shell as a function of time. At the current 
radius of the bubble in G298.4-0.3, r ^ 55pc, the radiation 
has deposited about twice the momentum supplied by the HII 
gas pressure. The stellar jets are not active at this time, but 
over the time they were active (corresponding to radii below 
^ 15pc) they deposited a momentum comparable to that of 
the radiation pressure (at those early times). In these mod- 
els, a combination of proto-stellar jets and radiation pressure 
disrupts the natal cluster, while a combination of gas and ra- 
diation pressure disrupts the GMC. The shell velocity at late 
times is of the order of the turbulent velocity seen in the ISM 
of the Galaxy (upper right panel), demonstrating that, even in 
the absence of supernovae, massive star formation can gener- 
ate turbulent motions on large (50pc or larger) scales compa- 
rable to those observed. 

4.2. The Starburst M82 

M8 2 is one of the nearest (D = 3.6 Mpc; iFreedman et alj 
119941) starburst galaxies, with an i nfrared luminosity Lir = 
5.8 X 1O'°L0 jSanders et all 12003b . The galaxy is small 
compared to the Milky Way, with a circular velocity Vc ~ 
llOkms"' (Yo ung & Scovilldll984 , and a CO i nferre d gas 
mass 2 x K^Mq inside r = 350pc dWeiB et all llOOl ) (ad- 
justed to our assumed distance), yielding a gas surface density 
Y,g « O.lgcm"^ and a gas fraction k, 0.2. The metallicity 
is 1.2-2.0 times solar (Smith et al. 2006). 

The radius and mass of the most massive star clusters in 
M82 are well established; there are about 200 clusters with 
M > lO^M© (Melo et al. 2005) and about - 20 well stud- 
ied super star clusters (Md > IO^Mq). With one arcsec- 



ond corresponding to a spatial scale of 17.5 pc, a number of 
super star clusters are resolved by HST ( Smith & Gallaghe^ 
l200ltlMcCradv et al.ll2003l:lMcCradv & Graham.2007) . typ- 
ical half Ught proj ected radii for these m assive objects are 
- 0.08" or 1.4pc. ' McCradv et all (12003 ') Hst 20 such clus- 
ters. The total mass of the 15 clusters for which they measure 
viral masses is ~ 1.4 x lO^M©. Their largest cluster, 'L', is a 
monster, with a mass of 4 x IQ^Mq and a half-light radius of 
1.5pc; more typical masses are ^ 5 x 10^M(D. A rough fit of 
the form (|4| gives « 1.9 (iMcCradv & Grahamlf2007i) . 

The masses of the GMCs in M82 are also known; the distri- 
bution is well fitted by equation with a n» 1.5 ±().l, an d 
a maximum mass of Mgmc « 3 x 10''Mq keto et alj |2005'). 
The Toomre mass is « 7 x lO^M©. Both are comparable 
to the mass of the two la rgest super star clusters given by 
iMcCradv & GrahamI (l2007i) . Either there were more massive 
GMCs in M82 in the past, or egmc ~ 1 for the GMCs out of 
which these two clusters formed. 

Tables 1 & 2 summarize our assumed galaxy, GMC, and 
star cluster properties in M82. These are all motivated by, 
and reasonably consistent with, the observations summarized 
above. Our results for the disruption of GMCs are summa- 
rized in Figure |2l The top two panels show the shell radius 
as a function of time and the velocity as a function of radius. 
The main sequence lifetime of a 1 20Mq star (the dashed ver- 
tical line) is comparable to the disk dynamical time (R^i / vv, 
the solid vertical line). The GMC is disrupted (reaches the 
Hill radius) about one disk dynamical time after the cluster 
forms. The velocity of the shell reaches higher values than 
those found in our Milky Way model because of the much 
larger cluster masses in M82, combined with the fact that the 
star cluster radii in the two galaxies are nearly the same. Ini- 
tially, the shell velocity is comparable to the escape velocity 
from the cluster The velocity begins to slow once the swept 
up mass is similar to the mass in the cluster (at r ~ 4pc). 

What is responsible for disrupting the GMC in our M82 
model? The lower panels of Figure |2] show that the cluster 
gas on small scales is expelled by a combination of proto- 
stellar jets and radiation pressure, while the overlying GMC 
is disrupted primarily by radiation pressure. The contribution 
from the gas pressure of the HII region is negligible over most 
of the evolution: the lower right panel shows that the H II gas 
pressure contribution to the momentum is ^ 10% of that con- 
tributed by radiation when the shell radius reaches the original 
size of the GMC. 

4.3. The Clump Galaxy Q2346-BX 482 

Q2346-BX 482 is a redshift z = 2.26 disk galaxy with a 
disk radius of Rd « 7kpc and a gas mass, as estimated from 
inverting the Schmidt Law in Kennicutt (1998), of Mg w 
3 X 1O'°M0 dGenzel et al.ll2008l) . We interpret the clumps in 
rapidly star forming redshift z ~ 2 galaxies as Toomre mass 
GMCs, with radii Rgmc ~ Ikpc, and we model the giant 
clump in BX 482, as one of the most extreme examples of 
this phenomenon. With turbulent and circular velocities of 
vj ~ 55 km s~' and vv « 235 km s~', respectively, we infer a 
disk scale height of // w 1.6kpc, and using 0g = 3 a GMC 
size of order 500 pc. For comparison, the observed clumps 
are modestly larger than this, around 1 kpc. 

The mass of the central star cluster is set so that the lumi- 
nosity matches that observed, roughly L w 4 x lO^L© (Gen- 
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Fig. 2. — Shell radius as a function of time and velocity, forces, and shell momentum as a function of radius in our model for the disruption of a GMC by a 
massive star cluster in the starburst M82. Upper left: Bubble radius as a function of time. The dashed line is when the star cluster luminosity drops to 1/3 of its 
initial value, while the dotted vertical Hne marks the lifetime of a 120Mq star. The shell is disrupted by ambient differential rotation and tidal forces in the galaxy 
at the Hill radius (solid line); this occurs roughly when the first supemovae explode. Upper right: The velocity of the swept-up shell in the M82 model as a 
function of shell radius. Lower left: The forces in our M82 model, plotted against the radius of the swept up shell. Line styles are the same as in Figure[T] Lower 
right: The momentum of the swept-up shell in the M82 model, with the same line styles as in Figure[T| The bulk of the momentum is supplied by radiation; the 
contribution from gas pressure is negligible. The early contribution of the protostellar jets is important in the disruption of the natal cluster gas. 



zel, private communication). We take the radius of the star 
cluster to be ~ 7pc (see ^2.41 ). In reahty, there will likely be 
many star clusters, with a distribution given by equationlH and 
with a spread in age of several to ten Myrs. We assume solar 
metallicity, consistent with the observations. Finally, the gas 
mass of the clump is not known, but we assume it is roughly 
the Toomre mass, ^ IO^Mq. This is consistent with the mass 
of ionized gas for the observed luminosity of the clump, at the 
observed size Rgmc ~lkpc, and for a stellar population less 
than 4 Myrs old (eqs.lAT0l& lAT2l) . 

Figure [3] shows the results of our model for the giant GMC 
in Q2346-BX 482. The right panel shows that radiation pres- 
sure dominates the evolution of the GMC at nearly all times. 
The GMC is disrupted in about 15 Myrs, half the disk dynam- 
ical time scale. The shell velocity, shown in the middle panel. 



is ^ 30- 80km s~ when the radius is 0.5 - 1 kpc, in reason- 
able agreement with the observed velocity dispersion of the 
galaxy. Note that the decrease in velocity at late times (due 
to the decrease in radiation pressure seen in the right panel) is 
probably not that physical given that the shell is in the process 
of being sheared apart by the rotation of the galaxy. 

Our conclusion that radiation pressure is disrupting the 
massive clump in BX482 is directly supported by observa- 
tions, independent of the specific assumptions in our model: 
the self-gravity of the clump is 

GMM 34/ M \Vlkpc\' 

^..v = ^ = 3xlO (^^j (^^j dynes, 

(10) 

where we have scaled M to the Toomre mass. This can be 
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compared directly with the radiation pressure force. 



34 



F,ad=-=5xlO 
c 



The clump should thus be expanding. 



dynes. 



(11) 



4.4. The Ultra-Luminous Infrared Galaxy Arp 220 

Arp 220, at 77 Mpc, is the proto-typical ULIRG in the 
local universe. The gas mass of each of the w lOOpc 
star-forming disks in Arp 220 is 1O^M0, the circular velocity 
Vc ~ 300km s~', v j » 80km s~' , and the disk scale height H = 
(vt I Vc)'"<i ~ 23 pc (iDownes & Solomonll998tlSakamoto et al.l 
[1999'). The mean surface density is Ej « 7 gem , about 100 
times larger than that of M82 and more than a thousand times 
higher than in the Milky Way. We estimate that Arp 220 has 
GMC masses of « 5 x 1O^M0, Rgmc ~ 5pc and a turbulent 
velocity in each GMC of ~ 170km s~', about twice that mea- 
sured for the disk as a whole. Although the metallicity is un- 
certain, we take a fiducial metallicity of 3 times solar; this in- 
crease in metallicity is important because it increases the dust 
optical depth and hence the overall importance of radiation 
pressure. 

IWilson et al.l (l2006h found ^ 40 young superstar clusters in 
and around Arp 220; they estimate masses for about a dozen, 
with a number having 2-4 x 1O''M0; the largest has 

Mci « IO^Mq. Given the huge extinction toward the twin 
disks, this is likely to be a rather conservative lower limit on 
the mass of the most massive cluster in the system. The clus- 
ters are unresolved in the HST images (which have a resolu- 
tion of order 15pc at the distance of Arp 220), except possi- 
bly thei r brightest cluster, with a half light radius r^\ « 20 pc. 
I Wilson et al. ( 2006 ) do not obtain either a velocity dispersion 
or a half light radius for their clusters, so they cannot calcu- 
late a dynam ical mass. Rather, they use a Salpeter IMF and 
iBruzual & C harlot (1993) stellar synthesis models combined 
with their photometry. 

We find that for Md = 1 .4 x IO^Mq (L = 3 x lO'^L©), even 
a Toomre mass GMC (4 x WMq) in Arp 220 would be dis- 
rupted (see Fig. |4]i. The disruption of the GMC occurs on 
the dynamical time of the disk, well before any supernovae 
explode in the GMC's central star clusters. 

Our estimated GMC mass in Arp 220 (5 x WMq) is a fac- 
tor ten higher than the largest GMC mass seen in M82; this 
is a result of the much larger surface density in Arp 220 com- 
pared to M82. In contrast, the star cluster masses found so 
far in Arp 220 are only a factor 2-3 times larger than the 
masses of the largest clust ers observed in M82, the la tter be- 
ing around 2-4 x \&Mq dMcCradv & Grahaml20()7l) . Given 
that our predicted GMC star formation efficiency is not that 
different in the two cases (Table 2), we expect that more mas- 
sive clusters are lurking in Arp 220. 

The right panel of Figure |4] shows the forces acting on the 
shell of swept up mass in our model of Arp 220. As in Figures 
[T|and|2] the force due to proto-stellar jets is initially similar to 
that due to radiation pressure. This situation lasts only while 
the shell accelerates from the initial clump radius of about 
4pc, until the shell reaches a little less than 6pc. For the rest 
of the evolution radiation pressure provides the only signif- 
icant outward force. The outward force supplied by ionized 
gas is completely negligible; the short dash line in the Figure 



is the gas pressure multiplied by 100. Both the hot gas and 
cosmic rays produces by shocked stellar winds are dynam- 
ically unimportant, even though we have assumed complete 
trapping of the shocked stellar wind material (an assumption 
that fails in the Milky Way; Harper-Clark & Murray 2009). 

The middle panel of Figure |4] shows that radiation pressure 
will stir the ISM of Arp 220 to ~ 50km s~', somewhat less 
than the escape velocity from the cluster and similar to the 
velocity dispersion seen in CO observations. 

5. DISCUSSION 

5.1. The Importance of Radiation Pressure 

Using four examples that cover conditions ranging from 
Milky Way-like spirals to the densest starbursts (see Tables 
1 & 2), we have explored the physical processes that can dis- 
rupt giant molecular clouds (GMCs), one of the basic build- 
ing blocks of star formation. We find that radiation pres- 
sure produced by the absorption and scattering of starlight by 
dust grains can contribute significantly to disrupting GMCs in 
nearly all types of galaxies. By contrast, protostellar jets are 
important only at early times during GMC disruption while 
the thermal gas pressure in HII regions is important for GMC 
dispersal in spiral galaxies like the Milky Way, but not in more 
luminous starbursts. For the Milky Way and M82, where the 
observations are particularly detailed, our results demonstrate 
that observed massive star clusters have precisely the lumi- 
nosities and structural properties required to disrupt Toomre- 
mass GMCs via radiation pressure. 

The results presented here support Thompson et al. (2005)'s 
model of radiation pressure supported star-forming galaxies. 
In that paper, we focused on the large-scale properties of star 
formation in galaxies and the fueling of massive black holes in 
galactic nuclei. Here we have extended that model by taking 
into account the fact that star formation is not smooth and 
homogeneous; rather, most stars form in massive star clusters 
inside massive GMCs (©. Our conclusion that Toomre-mass 
GMCs can be disrupted by radiation pressure is qualitatively 
and quantitatively similar to Thompson et al.'s conclusion that 
radiation pressure can regulate star formation in galactic disks 
to have Toomre's Q~\. 

It is useful to consider simple scaling arguments in order to 
understand why, over the range of surface densities probed by 
the observed Schmidt Law (10~^ g cm"^ ^^g% 10 g cm"^), 
radiation pressure is the most viable mechanism for GMC dis- 
ruption. To rough approximation, the self-gravity of the gas 
in a GMC is 
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which varies by a more than a factor of ^ 10* from normal 
galaxies to starbursts. 

We can compare this force directly to the radiation pressure 
force. In the optically thin hmit^ 
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By this we mean that the GMC is optically-thick to the UV, but optically- 
thin to the re-radiated FIR emission. 



10 



0.01 



10 



t (Myrs) 




Fig. 3. — Left: The radius of the shell as a function of time in the model for the giant clump in the z = 2.26 galaxy Q2346-BX 482. The dotted, dashed, and 
solid lines mark when the first cluster SNe explode, the central cluster luminosity drops to 1 /3 of its initial value, and t reaches the disk's dynamical timescale, 
respectively. Middle: The velocity of the swept-up shell. For the vast majority of the evolution, the shell velocity is comparable to the observed velocity dispersion 
of the gas, ~ 55 km s"' jGenzel et alJ2008i) . Right: The forces in our model for the giant clump in BX482, with line styles as in Figure[T] 
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Fig. 4. — Left: Bubble radius as a function of time in a model for the disruption of a GMC by a star cluster in the ULIRG Arp 220. The dotted, dashed, and 
solid lines mark when the first cluster SNe explode, the central cluster luminosity drops to 1 /3 of its initial value, and 1 reaches the disk's dynamical timescale, 
respectively. Note that in Arp 220, the disk and GMC dynamical times are short compared to the main sequence lifetime of massive stars, unlike in our other 
models (see also Table 1). Middle: The velocity of the swept-up shell as a function of shell radius. Right: The forces as a function of the radius of the swept 
up shell. The upper most solid line is the total outward force Fom, consisting of five components: the force exerted by protostellar jets (solid line; fjei). the force 
exerted by HII gas pressure multiplied hy 100 (short dash line; Fhh) the force exerted by radiation pressure on dust grains (long dash line; /Vad). the force exerted 
by shocked stellar winds (Fjo,; soUd), and the force exerted by cosmic-rays produced in stellar wind shocks (F„ ; dot-dashed). Radiation pressure dominates the 
outward force at nearly all times. The dotted line is the total inward force Fi„, dominated by the self-gravity of the shell. 



where vj/ is the light-to-mass ratio in cgs units. Thus, 

4>IGcMt /Mt\ 

^ /sooocgsN (om\ / ^ \ 

I ^ VecMcy Uxl0-3gcm-2;'^'^^ 

where we have scaled to values appropriate to the Galaxy. We 
see that if egmc increases with gas surface density, as our 
calculations indicate (e.g.. Table 2), then radiation pressure 
provides a plausible mechanism for GMC disruption in both 
spiral galaxies like the Milky Way and somewhat denser and 
more luminous galaxies. 

For galaxies with sufficiently large surface densities, > 
0.5 g cm"^, GMCs will be opaque to the emission by dust in 
the far-infrared. This increases the radiation pressure force so 



that (in the thin shell approximation used here) 

L 

Frad = Ti-ad - OC Sg , (15) 

where T,-ad = KFiRSsh/2, kfir is the Rosseland mean opacity of 
the GMC in the FIR and Sjh is the surface density of material 
in the shell. Comparing the optically thick radiation pressure 
force with that due to gravity, we find that 

fih 4ttGc /Mr\ 

^ ^ /3000cgs\ / 30cm^g-' \ / 0.25 \ ^^^^ 

where we have scaled to a relatively high value for the 
Rosseland-mean dust opacity (see below). Note that the ra- 
tio in equation ( fT6] l does not explicitly depend on stellar/gas 
mass, because both and Fj-ad are oc M^. Using the scalings 
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in Appendix IaI it is easy to see that no other stellar feedback 
process has this property. Indeed, most of the previously sug- 
gested support mechanisms scale as F oc with/3< 1, VIZ, 
HIl gas pressure, stellar winds, and pressures associated with 
shocked stellar winds. For this reason, although many feed- 
back mechanisms are competitive with gravity in normal spi- 
rals, the self-gravity of the disk quickly overwhelms the forces 
due to stellar feedback in starburst galaxies. In contrast to 
these other feedback processes, radiation pressure in optically 
thick gas scales as Fi-^^ oc M^Mg, so that it is at least in prin- 
ciple possible that radiation pressure can disrupt GMCs even 
in the densest, most gas-rich environments (e.g., ULIRGs and 
z^2 galaxies). Radiation pressure is, to our knowledge, the 
only stellar feedback process that has this property. 

Equation ( fTSI l shows that the efficiency of star formation 
in a GMC at very high densities is sensitive to the metal- 
licity and dust composition, which influence the FIR opac- 
ity KpiR, and to the stellar IMF, which determines the light 
to mass ratio of the stellar population 5*. In very dense en- 
vironments there are some reasons for su specting that the 
IMF may be top heavy (e.g.. ' Murravll200 9*), as appears to be 
the case in regions of massive star formation more generally 
(iKrumholz et al.ll2007l) . If this is indeed the case, it would 
increase ^ and thus decrease the ecivic required for GMC dis- 
ruption.^ For a relatively normal IMF, however, the star for- 
mation efficiency in GMCs must be appreciable at high den- 
sities, with ecMC ~ 0.25 or perhaps even larger. 

5.2. The FIR Optical Depth 

A key part of our argument for the importance of radiation 
pressure is the fact that star clusters have very high surface 
densities, which can trap the FIR radiation produced by dust, 
thus enhancing the radiation force by a factor of ^ Trad in the 
optically-thick limit (eq. |[T6l ). Figure |5] shows the Rosse- 
land mean dust optical depth (T,ad) through the shell as a func- 
tion of radius in the M82 and Arp 220 models (compare with 
Figs.m&IDi. In the M82 model, r^ad > 1 for r < 20pc, while 
for Arp 220 Trad 3> 1 at all radii. In our Milky Way models, 
by contrast, we find that the GMC is essentially always in the 
optically-thin limit, i.e., opaque to the UV but not to the FIR. 
Figure|5]also shows the effective optical depth in both the M82 
and Arp 220 models (dashed lines), which we define as 



Teff : 



(Lt/c) 



(17) 



The effective optical depth quantifies the enhanced coupling 
of photons emitted by a cluster in the center of a GMC com- 
pared to photons originating at the mid-plane of a uniform 
density disk. In M82 the effective optical depth at the end 
of the bubble evolution is about equal to the mean optical 
depth at the mid-plane. By contrast, in our Arp 220 mod- 
els, the momentum deposited per photon in the bubble shell 
is a factor of ~ 3 larger than would deposited by a photon 
traversing a uniform density disk. This shows that the effect 
of the radiation pressure may be three times larger than cal- 
cu lated using the mean m id-plane optical depth, as was done 
in [Thompson et al.l (l2005h . This effect increases the impor- 
tance of radiation pressure in the densest galaxies, where it 

' Note that as IMF becomes arbitrarily top-heavy ^ — > 4ttGc/kt, 
where kt is the Thomson opacity. This sets a minimum on the ratio 
Fsh/Fy^a for any stellar population: -Fsh/ZVadLin ^ (kt /kfir) Sgmc ~ 
10-2(30cm2g-VKFiR)eo'^^c- 



is needed most, effectively decreasing the normalization in 
equation ( fT6] l. Note also that because most of the star for- 
mation — and thus radiation — occurs in a few massive star 
clusters in the most massive GMCs (§|2]i, there is very little 
"cancellation" due to different radiation sources drivi ng the 
gas in different directions (as was suggested by Socrat es et al.l 
120061) : the distribution of radiation sources in galaxies is not 
well-approximated as infinite and homogeneous (§|2]l. 

Given the turbulent and clumpy nature of the ISM in GMCs, 
one may question whether or not the photon coupling ef- 
ficiency is as large as Trad or Tgff, since these expressions 
assume uniform shells of matter We have, after all, used 
the argument that GMCs are porous to argue that hot gas 
from shocked stellar winds escapes rapidly from GMCs in the 
Galaxy. The optical depth r^d is measured from the center 
of the protocluster outward and is proportional to the column 
density of overlying gas; the latter has been observationally 
measured by a nu mber of authors (e.g., Goodman et al.ll2009i: 
IWong et al.ll2008t) and is consistent with a log-normal distri- 
bution. Numerical si mulations also find log-normal surface 
density distributions (Ostriker et al."2001).'" In the notation 
of Ostriker et al. (2001), Goodman et al. (2009) find a disper- 
sion in the logarithm of the column density of 0.11 < a < 
0.22, corresponding to a range of mean logarithmic column 
densities 0.01 < fj, < 0.05. This agre es well with the high 
Mach number turbulence simulations of lOstriker et al.l(l2001l) . 
For /i = 0.05, 99% of sight lines have r/f > 0.5, where f is the 
(angular) mean of the optical depth. There are thus very few 
optically thin sight lines until the mean optical depth is itself 
of order unity. This conclusion is based on observations that 
do not probe, and simulations that do not include, some of the 
physics relevant to our models (e.g., the Rayleigh-Taylor and 
photon-bubble instabilities; Blaes & Socrates 2003; Turner et 
al. 2007; Thompson 2008; Krumholz et al. 2009). Nonethe- 
less, these results provide some support for our argument that 
the dense ISM will be opaque to FIR radiation, thus signifi- 
cantly enhancing the radiation pressure force in the optically 
thick limit. In the absence of very large changes in the IMF 
(^P) or the average dust-to-gas ratio (kfir), this enhancement 
is required for radiation pressure to disrupt GMCs in dense 
starbursts. For Milky Way-like galaxies, our conclusions are 
not sensitive to the radiation force in the optically thick limit, 
since the GMCs/star clusters are not optically-thick except on 
very small scales. 



5.3. The Origin of Large-Scale Turbulent Motions 

In all of our calculations, the shell velocity at r w Z/^/ is sim- 
ilar to the turbulent velocity required to maintain Toomre's 
2 ~ 1 in the ambient disk, thus staving off self-gravity on 
the largest scales. This is not a coincidence: disrupting the 
GMC — which we have explicitly taken to be the Toomre 
mass — requires a force comparable to that needed to stir up 
the gas in the disk proper to 2 ^ 1 • We have interpreted the 
shell velocity at large radii as a turbulent velocity because the 
disruption of the shell by shear, tidal forces, and instabilities 
(e.g., Rayleigh-Taylor) will inevitably convert this relatively 
ordered kinetic energy into random fluid motions. More so- 
phisticated models are clearly called for, but it is also clear 

The observations measure Traj along the line of sight from the Earth 
through the cloud rather than from the center of the cloud outward, but the 
two surface density distributions should not differ dramatically. 
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Fig. 5. — The Rosseland mean optical depth (1 +rrad) (soUd lines) and the effective optical depth Teff (eq. 117) (dashed lines) as a function of radius in the M82 
(left) and Arp 220 (right) models, respectively. Note that Teff is a factor ~ 3 larger than the mean optical depth of the disk in Arp 220 (the horizontal dotted line), 
enhancing the magnitude of the radiation pressure force. 



that such models must include the effects of radiation pres- 
sure, particularly in very optically thick systems like Arp 220, 
but also in models of Milky Way-like galaxies. 

5.4. The Global Star Formation Efficiency in Galaxies 

A key question we are unable to fully address is how the 
efficiency of star formation in a GMC (egmc) relates to the 
global star formation rate in the galaxy as a whole, as encap- 
sulated by, e.g., the parameter 77 in the Kennicutt Law (eq. IT]) 
or the observed Schmidt Law. We briefly discuss some of the 
relevant issues here, but leave a more detailed analysis of this 
problem to future work. 

In our models of GMCs in Milky- Way like spirals, the time 
for a star cluster to disrupt the natal GMC is 1/3 to 1/2 the 
local dynamical time. When we add to this the time for the 
GMC to contract to its present size, the total time involved 
will be 2 longer than the local dynamical time of the disk. If 
the subsequent supernovae do not greatly prolong the process 
of reincorporating the bulk of the gas back into the disk, which 
we suspect is correct, the time averaged star formation rate per 
unit area will be 

= eGMcS/(2fdyn) = JyS/fdyn, (18) 

with 7/ given roughly hy rj ^ eGMc/2 ~ 0.02, in good agree- 
ment with the observations. 

The number of clusters capable of disrupting a Toomre 
mass GMC, in our feedback model, should be proportional to 
the number of Toomre masses, ~ (R/H)- ^ 700. The corre- 
sponding number of giant Hll regions observable at any time 
is {R/HftMs/(2tdyn) ~ 30, where fMS « 3.6 x 10* years is the 
lifetime of an early O star, which is also the time over which 
a star cluster will emit a large luminosity of ionizing photons. 



This estimate is consistent with the fact that half of the free- 
free emi ssion in the Galaxy is prod uced by about 17 giant HII 
regions dMurrav & Rahmanll2009l) . 

The GMC star formation efficiency is eoivic = 0.38 for 
our fiducial Arp 220 model (Table 2); more generally, it is 
eoMC ^ 0.25 in the optically thick limit for typical IMFs 
and dust to gas ratios (eq. |[T6l ). This suggests that 77 is 
also « 0.25, higher than implied by the Kennicutt relation 
(eq. IXI), although reasonably consistent with the conclusions 
of Bouche et al. (2007). However, this estimate assumes that 
after the dispersal of the GMC, gas falls back into the disk 
and forms a new GMC in a single dynamical time; as the left 
panel of Figure |4] shows, the timescale for cluster disruption 
is in fact ^ 1/10 of the timescale for the cluster luminosity 
to drop significantly. Thus, the star cluster luminosity can 
drive motions in the remaining gas ^ vj ^ {H/R)vc, suffi- 
cient for hydrostatic equilibrium over many ('^ 10) dynamical 
timescales. Because the photon diffusion timescale is rapid 
compared to the dynamical timescale, the medium cannot be 
supported stably (Thompson et al. 2005; Thompson 2008). 
Hydrostatic balance will only be maintained in a statistical 
sense within a volume ^ 4H^ of the star cluster. After the 
cluster's luminosity decreases on a timescale ^ Jms the gas 
will recollapse to form a new GMC and the process will re- 
peat until gas exhaustion. We suspect that this may reduce r/ 
by a factor of ~ fdyn/fMS, but clearly more work is needed. 
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APPENDIX 

FORCES INCLUDED IN THE MODELS 
In this section we describe the forces acting on gas surrounding massive star clusters embedded in giant molecular clouds. 

Inward Forces 
Gravity 

We assume that the gravitational force acting on the bubble shell consists of three components, Fgray = ^'stars + ^sheii+^disk- The 
central star cluster exerts a force on the bubble shell given by 

^;tars= 5 , (Al) 



the shell self-gravity is 

while the mass in the galactic disk exerts a force 



^sheU = -^, (A2) 



2r2 ' 

i^.s. = -^f. (A3) 

The last force is that exerted on the part of the shell rising vertically away from the disk; gas in the plane of the disk will not feel 
this force, but we ignore this complication, just as we ignore Coriolis forces. 

Ram Pressure 

As the shell expands into the ISM of the galaxy it will sweep up gas. This swept up gas exerts a ram pressure on the shell 
giving a force 

^;ain= 1— V^ (A4) 

ar 

where dM/dr is given by the appropriate derivative of equation dTji or ([8]l. 



Turbulent Pressure 



The pressure in the interstellar medium and in the GMC provides an inward force on the shell. We refer to these pressures as 
turbulent pressure, although there may be other components such as magnetic or cosmic ray pressure. We approximate them as 



^lui-b = -47rr PisM- 



r)2 

"gmc 



1- 



^GMC 



(A5) 



Here we have assumed that the GMC has a Larson-like mass distribution. A similar expression holds for clouds with p{r) ~ l/r^. 
The pressure of the ISM is given by 

PiSM ~ TTG^d- (A6) 



Outward Forces 
Gas Pressure Forces 

The large luminosity Q (number per second) of ionizing photons produced by the massive stars in a young star cluster will 
photoionize and heat gas in the vicinity of the cluster, raising the gas pressure above that in the neutral gas. This hot gas will 
exert an outward force on the bubble wall, with a magnitude 

FHii = 47rr2pHii (A7) 

The gas pressure is taken to be 

Phii = nkT (A8) 

where T = 8000 K ([Gardner et al.lll970h and 




(A9) 



The volume V = Anr^ /3 and the recombination coefficient arec « 4 x 10 In models for clusters in the Milky Way, the shell 
velocity is subsonic, so the gas pressure is roughly constant throughout the bubble. If there is any hot gas or cosmic rays due to. 
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for example, shocked stellar winds, they will also be roughly isobaric, and in pressure equilibrium with the HII gas. If the hot 
gas or cosmic ray pressure is in excess of the estimate given here, the HII gas will be confined to a fraction of the bubble volume. 
However, obse rvations of Carina and rnassive Milky Way and LMC clusters suggest that the pressure is well approximated by 
dAST i and (iHaroer-Clark & Murravll2009h . 

For large enough stellar clusters, those with at least one 35Mq star, Q is proportional to L, 

^=e« 8x10-" erg, (AlO) 

or about 3.6 Rydbergs (IRyd w 2.2 x 10"" is the energy required to ionize hydrogen). The HII gas pressure inside the bubble is 
then 

'■™='.ox >«-"'(t^) {^) (A") 

assuming a unit filling factor for the HII gas, i.e., ignoring the hot shocked winds. 
It is helpful to have an estimate for n in terms of the luminosity. The number density n(L, r) is given by 



"^^'^^=\/4;^- ^^^^^ 

If the bubble is breached, so that HII gas leaks out, the pressure in the bubble will be lower than the estimate given here, but the 
force experienced by the bubble wall will actually be larger The reason is that fewer ionizing photons are absorbed in the bubble 
cavity, leaving more to heat the gas on the interior of the bubble wall. The heated gas escapes away from the wall into a partial 
vacuum, exerting a force 

F=McHii (A13) 

on the wall. The mass loss rate is given by Anr^nipncii. In our numerical work we have assumed that the HII gas leaks out of the 
bubble, since the rapidly expanding ionized gas at the ionization front will exert a larger force, and we want an upper limit on the 
efficacy of HII gas pressure. 

In models for ULIRGs like Arp 220, the bubble expansion velocity is larger than the sound speed of 8000 K gas, so the HII 
pressure may be twice that in the subsonic case even if no gas escapes from the bubble. However, the pressure associated with 
HII gas is negligible in ULIRGs, as we show in the main text. 

Forces Associated With Shocked Stellar Winds 

Hot gas from shocked stellar winds is ofte n thought to be importan t in the formation of bubbles around massive stars or star 
clusters. We argue here that it is not; see also lHarper-Clark & Mu rray (2009). We further argue that cosmic ray pressure cannot 
be important in Milky Way bubbles; if they were, the bubbles would expand more rapidly than is observed. We show in this paper 
that neither hot gas nor cosmic ray pressure is not relevant in ULIRGs. 

The force due to hot, shocked stellar wind gas: X-Ray Constraints 

O stars emit high velocity ma s sive winds; in clusters these win ds are seen to shock and emit diffuse x-rays at ^ KeV energies 
(ISewardetalJll979t IOeylll996t ISmith. Struck. & Nowi3l2005h . If t he stellar winds ar e confined to the bubble interior, the 
associated pressure can be far larger than the ram pressure of the wind jCastor et alj|1975h . The force is given by 

Fh=ATTr^Ph, (A14) 
where Pi, = 2£'/,/(3y). The energy equation for the hot gas is 

dEh 

= U.-4TTr^PhV,k-AnlV; (A15) 

recall that V is the volume of a bubble of radius r and A is the cooling function. The wind luminosity L„ = (1 /2)77(Voc,/c)Lboh 
(77 « 0.5 is the fraction of stellar luminosity scattered in the wind) and Lboi is the bolometric luminosity of the cluster 

We follow ICastor et all (Il975h to find the number density «/, of hot gas, i.e., we assume that heat conduction (at the Spitzer rate 

5/2 

CspTjj ) drives a mass flow into the bubble interior (possibly supplemented from cold gas clouds in the interior of the bubble) at 
a rate 

dMh _ l6m7ipCspT^^^r 
dt ~ 25k ' 

Given E/, and the radius of the bubble we find the hot gas pressure; from equation (IAI6I 1 we find «/„ and hence the temperature. 
With «/, and 7},, we can find the cooling rate (the third term on the right hand side of the energy equation (IA15b . 

The solar-mass stars in clusters also emit x-rays; in low resolution (non-Chandra) observations, these stars appear as diffuse 
emission, so care must be taken to distinguish the two. Assuming that one can extract the stellar emission, the diffuse x-ray flux 
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can be used to constrain the pressure of any hot gas component. iHarper-Clark & Murray! (l2009h use this to show that for clusters 
in the Milky Way and the LMC, the x-ray gas is in approximate pressure equilibrium with the HIl gas, and that the Hll gas has 
a filling factor near unity (greater than ~ 0.1). The implication is that the shocked winds either radiate away the bulk of their 
energy, or else leak out of the bubble wall. In either ca se, they do not play any role in the dynamics of the bubble, a conclusion 
also reached bv lMcKee. Van Buren. & Lazaref5(ll984l) . 

IHarper-Clark & Murravl (120091) show that observations of ultracompact Hll regions dRauw et aLllIOOl iTsuiimoto et al.ll2006h 
reveal a similar story. 

If the bubbles that form around star clusters in ULlRGs are also leaky, the shocked wind pressure would be negligible, as in 
the Milky Way. However, even if we assume that the winds are perfectly trapped, the high ambient pressures in ULlRGs ensure 
that the shocked gas is not dynamically important. 

To see this, note that we expect a total number of young (ionizing) clusters in a ULIRG to be 




(A17) 



If there are ^ 20 clusters powering Arp 220, each must have Lci ~ 5 x lO'^L©, and a mass « 2 x IQ^Mq assuming a normal 
(iMuench et al. 2002) IMF; this cluster luminosity is also that needed to disrupt a GMC in Arp 220, as we showed in ^4.41 The 
wind luminosity of a single such cluster i s L„ Kil x 10'^' erg s~'. We note that Lci ~ 5 x lO'" is about 3-4 times the luminosities 
of known young star clusters in Arp 220 (I Wilson et al.ll2006b . 

From Figure|4]the bubble disrupts the GMC (^?gmc ~ lOpc) after a time f ~ 5 x 10'' s. Assuming no losses, the wind energy 
accumulated in the bubble is Eh = 2.5 x lO^^'erg. The pressure is 2Eh/3V = 10"^ dynescm"^, about equal to the ambient pressure 
in the disk but well below the radiation pressure, tF /c « 10~^dynescm"^. 

The force we have estimated from the hot gas is 1.6 x 10^^ dynes; in the right panel of Figure|4]the force from the hot gas is a 
factor of ~ 5 below this estimate. The reason is that the hot gas suffers radiative losses, so that we have overestimated £/,. We 
note that in our numerical work we have assumed that the hot gas fills the volume of the bubble; if it does not, the cooling rate 
will be higher than we have calculated. 

The force due to wind generated cosmic rays 

In addition to producing hot gas, wind shocks will generate cosmic rays. Perhaps a third of the wind energy may be expected 
to be deposited into cosmic rays. If, as observed in the Milky Way, the hot gas is advected out of the shell, it will probably take 
the cosmic rays with it, so that both the hot gas and the cosmic rays are not dynamically relevant. 

We note that the pressure due to cosmic rays, if perfectly confined, would be 1/6 that of similarly confined shocked winds. 
As noted in the introduction, the dynamics of bubbles in the Milky Way and the LMC are not consistent with the high pressures 
associated with trapped hot gas; it follows that cosmic rays are probably also not confined to the interiors of such bubbles. 

If the shocked winds in ULlRGs also flow through the swept up shell, cosmic rays are unlikely to be dynamically important. 
Even if the cosmic rays are confined, the cosmic ray pressure is still negligible. The cluster bolometric luminosity is 2 x 
lO^'^erg s~', so the wind luminosity is 5x10^' erg s"'. Assuming 1 /3 of this is deposited in cosmic rays, the cosmic ray luminosity 
is Lcr ~ 1 .6 X lO"'' erg s"' . The cosmic ray pressure when the bubble radius is lOpc is 

P^^ = ~ 7 X 10-" ( ' dynescm-2. (A18) 

This is about a half of the mean dynamical pressure in the galaxy, but only several percent of the dynamical pressure in the 
bubble. This is still likely to be an overestimate of the cosmic ray pressure, as cosmic rays will quickly diffuse out of the bubble. 
The cosmic ray mean free path in the Milky Way is 0.1 pc < < Ipc ( Schli ckeiser,.200Z) ; scahng to the lower value, the time 
for a cosmic ray to diffuse out of a GMC in Arp 220 is 

f„=%c^2xl0"f^')s, (A19) 

where = \tc/3 is the cosmic ray diffusion coefficient. This diffusion time is ~ 10-20 times shorter than the bubble dynamical 
time, so even if the mean free path for cosmic rays in Arp 220 is substantially smaller than in the Milky Way, cosmic rays will 
diffuse out of the bubbles in the ULIRG. Figure |4] shows that the cosmic ray pressure is a little less than one percent of the 
radiation pressure; the lower cosmic ray pressure is the result of taking Pci = L(^rtcr/V, which accounts for diffusive losses, instead 
of equation ( IA18b . Note that we have optimistically ignored cosmic ray losses due to inelastic proton-proton collisions, which 
cool the cosmic ray population with energies >GeV on a timescale fpp « 5 x 10''(10''cm-^/n)yr, where we have scaled to the 
average volumetric gas density of Arp 220. 

Protostellar Jets 

While stars are actively accreting in the protocluster, we assume that they will emit high velocity outflows or jets. Once the 
cluster disrupts, the accretion halts, and the jets turn off. We assume that this happens over a time given by 0ff times fff of the 
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cluster; in addition, we allow for an extra factor of two to account for the fact that the accretion disks will not vanish once the 
cluster gas is dissipated, but instead will accrete onto the star in a disk viscous time: 

fjet = 2(/)fffff . (A20) 
We take 0ff = 3. The jets deposit momentum into the surrounding gas at a rate 

^jet = ejetMaccVjet CXp-'/'J" (A21) 

The jets expel mass at a fraction ejet « 0.1 -0.3 of the mass accretion rate ( Matzner & McKee'2000'). at a velocity vjet that depends 
on the mass of the star. We estimate vjet by calculating the escape velocity from the surface of a star of mass m, using the radius 
of a star of that mass found from the main sequence radius as given by the Padova models; this is a slight over estimate of the 
escape velocity, since accreting stars are larger than main sequence stars, but we compensate for this by using the escape velocity 
rather than the (larger) terminal velocity of a wind or jet. We then average over the initial mass function to yield vjet. For a 
Muench et al. IMF we find vjet « 3 x 10^ cms"', but for top heavy IMFs this can increase to 10^ cms"'. We varied ejet over the 
ran ge 0.1-0.3 giving an ef fective velocity in the range 30 -90km s"', bracketing that found by, e.g.. ' Matzner & McKed (120001) 
and lNakamura & O ( l2007h . The results did not depend strongly on this parameter, but because the clusters we consider are so 
dense, we expect a substantial amount of momentum cancellation (from jets colliding nearly head on) we used low values for ejet 
in the runs presented here. 

Radiation Force 

The radiation force is given by 

f;ad = (l+Trad)-, (A22) 
C 

where Tiad is the Rosseland mean opacity through the shell. The gas in the shell has a temperature that is of order lOOK, so for 
most clusters in the Milky Way Tiad < 1. We calculate the optical depth as follows: we assume the shell has a thickness 1 /lO of 
its radius, and divide it into 100 pieces. We assume that the density i s constant through the shell. We calculate Teff from L and 
the radius r of the shell. Given T and p, we use the opacity tables of Se menov et al.l (l2003h to find the Rosseland mean opacity 
in the outermost slice of the shell. Subsequently we find the optical depth, then integrate the radiative transfer equations inward 
through the shell. The factor of unity in equation ( IA22l l accounts for the fact that the temperature of the radiation field is of order 
30,000K at the inside edge of the shell, before the photons have encountered any dust grains; upon striking a dust grain, the 
photons provide an impulse per unit time given by L/c. 



